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Abstract 

Recently f(T) modified teleparallel gravity where T is the torsion scalar 
has been proposed as the natural gravitational alternative for dark energy. 
We perform a detailed dynamical analysis of these models and find con- 
ditions for the cosmological viability of f(T) dark energy models as geo- 
metrical constraints on the derivatives of these models. We show that in 
the phase space exists two cosmologically viable trajectory which (i) The 
universe would start from an unstable radiation point, then pass a saddle 
standard matter point which is followed by accelerated expansion de sitter 
point, (ii) The universe starts from a saddle radiation epoch, then falls onto 
the stable matter era and the system can not evolve to the dark energy 
dominated epoch. Finally, for a number of f(T) dark energy models were 
proposed in the more literature, the viability conditions are investigated. 
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1 Introduction 



Recent cosmological observations indicate that our universe is in accelerated 
expansion. These observations are those which is obtained by SNe la [JJ, 
WMAP [2], SDSS [3] and X-ray g]. These observations also suggest that 
our universe is spatially flat, and consists of about 70% dark energy (DE) 
with negative pressure, 30% dust matter (cold dark matter plus baryons), 
and negligible radiation. In order to explain why the cosmic acceleration 
happens, many theories have been proposed. The simplest candidate of the 
dark energy is a tiny positive time-independent cosmological constant A, for 
which oj = — 1. However, it is difficult to understand why the cosmological 
constant is about 120 orders of magnitude smaller than its natural expecta- 
tion (the Planck energy density). This is the so-called cosmological constant 
problem. Another puzzle of the dark energy is the cosmological coincidence 
problem: why are we living in an epoch in which the dark energy density 
and the dust matter energy are comparable?. An alternative proposal for 
dark energy is the dynamical dark energy scenario. The dynamical nature 
of dark energy, at least in an effective level, can originate from various fields, 
such is a canonical scalar field (quintessence) [5], a phantom field, that is a 
scalar field with a negative sign of the kinetic term [6], or the combination 
of quintessence and phantom in a unified model named quintom [7]. In the 
other hand modified models of gravity provides the natural gravitational 
alternative for dark energy [8]. Moreover, modified gravity present natu- 
ral unification of the early time in action and late-time acceleration thanks 
to different role of gravitational terms relevant at small and at large cur- 
vature. Also modified gravity may naturally describe the transition from 
non-phantom phase to phantom one without necessity to introduce the ex- 
otic matter. Among these theories, scalar-tensor theories [9], f(R) gravity 
[TO] are studied extensively. 

Recently a modified gravitational theory, the f(T) gravity, has been 
proposed |13|. HH] to explain the acceleration of the cosmic expansion which 
attracts much attention. In the simple case f(T) = T, the f(T) theory 
can be directly reduced to the Teleparallel Equivalent of General Relativity 
(TEGR) which was first propounded by Einstein in 1928 [12] . Similar to 
the f(R) theories, f(T) theories deviate from Einstein gravity by a function 
f(T) in the Lagrangian, where T is the so-called torsion scalar. Types of 
these models have been proposed to explain the late-time acceleration of the 
cosmic expansion without including the exotic dark energy [13]-|27j. 

In this paper we will investigate the cosmologically viable conditions for 
general f(T) gravity. Wu et. al, [22] and Zhang et. al, [23] in a similar 
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way have investigated the dynamical behaviors of T + f(T) model. They 
obtained two critical points for radiation and matter dominated, respectively 
and a critical line for dark energy dominated. Then for a concrete power- 
law model, they obtained that dynamical system has a stable de sitter phase 
along with an unstable radiation dominated phase and an unstable matter 
dominated one. 

However, in our study for general f(T) model we obtain three critical 
continuous lines corresponding to de sitter points (O m = 0,Q r = 0,Qde = 
l,cj e ff = —1), scaling solution with matter (Cl m = l,£l r = 0,$1de = 
0,u e ff = 0) and scaling solution with radiation (f2 m = 0, Q r = 1, Qde = 
0,u) e ff = |), respectively, we find the general conditions for the cosmolog- 
ical viablility of these models as geometrical constraints on the derivatives 
their and obtain two kinds of cosmological viable trajectory for general f(T) 
model: (i) we have a saddle matter phase which is followed by a de sitter 
epoch along with an unstable radiation era (ii) the system will stay in the 
stable matter dominated phase and can not evolve to dark energy dominated 
which would start from a saddle radiation era. 

The paper is organized as follows, in the following section, we review 
f(T) theories and in order to study the dynamics of these models, we have 
to rewrite the field equations into a 4 dimensional autonomous system which 
once is a constant, using the method we construct the curve m(r) 
which m = -jmp- and r = — -jjjy , where a prime denotes derivative with 
respect to the torsion scalar T. In section [3l we drive the critical lines 
with their stabilities then in the (r, m) plane, we obtain conditions for the 
cosmological viable of f(T) dark energy models. Section [J] is based on the 
numerical analysis for a number of f(T) models to validity the analytical 
results presented in the section [3j Finally, in Section [5J we present our 
conclusions. 

2 f(T) dark energy models 

In this section firstly, we briefly review the Teleparallel gravity so-called f(T) 
gravity in the spatially flat FRW universe. Then we obtain the autonomous 
equations to study cosmological dynamics of f(T) gravity models. 

A. Field equations and definitions 

We start with the generic form of the action of Teleparallel gravity as 

5 = 4/ edx 4 (f(T) + L r + L m ), (1) 
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where k 2 = 87rQ, e = \J—g = det(e 1 ^) and T is the torsion scalar, also L r and 
L m are the Lagrangian density of the radiation and the matter, respectively. 

In the f(T) framework, in the spatially flat FRW universe e^, the veribein 
field is related to the metric as 

here p, u are the coordinate indices on the manifold while i,j are the coor- 
dinate indices for the tangent space of the manifold which all indices run 
over 0, 1,2,3, also r]ij=di&g(l, -1, -1, -1). 

Varying the action Eq.([T]) with respect to the veribein field leads to the 
equations [14] 



SrMT)fMT) + \e-%(eSn ~ MT) (3) 



1 k 
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Here a prime denotes differentiations with respect to the torsion scalar T. 
Also S^ v = e^Sp U and T^ u is the matter energy- momentum tensor. 
In the spatially flat FRW metric <7 Ml/ =diag(-l, a(t), a(t), a(t)), the set of 
field equations Eq.([3|) for i = = v reduce to |14j 



l2H 2 f'(T) + f(T) = 2k 2 (p r + Pm ). (4) 

and for i = 1 = v lead to 

A8H 2 H f" (T) - (12H 2 + AH)f(T) - f(T) = 2k 2 ( Pr +p m ). (5) 

and torsion scalar as a function of the Hubble parameter H = - 

T = —QH 2 . (6) 

Where a dot represents a derivative with respect to the cosmic time t also, 
(p r ,Pr) and (p 

miPm) the (total energy density, pressure) of the radia- 
tion and the matter inside the universe, respectively. 

Using Eqs.([U E]), one can obtain the modified Friedmann equations as 
follows 

3H 2 = k 2 (p r + p m + PT ), (7) 
2H + 3H 2 = -k 2 {p r + Pm +p T ), (8) 
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with 

PT = ^(2Tf'(T)-f(T)-T), (9) 

PT = -^H{-2Tf"{T)-f\T) + l)]-p T . (10) 

Here pt and pt are the torsion contributions to the total energy density and 
pressure, respectively. 

The continuity equations of these energy densities p r , p rn , px are: 

= pr + 4Hp r , (11) 
= p m + 3H Pm , (12) 
= pt + 3H(pt+Pt). (13) 

Prom Eqs.([9|), (10), we can define gravitationally induced form of dark 
energy density pt = Pde and pressure px = Pde- The equation of state 
parameter is defined as 

= PDE -8HTf"(T) + (2T - AH)f'(T) - f(T) + AH-T 
UDE Pde 2Tf'(T)-f(T)-T ' 1 j 

and for a f(T) (DE) dominated universe, one can obtain the effective equa- 
tion of state 

2H 2Tf"(T) + f'(T) - 1 
3^ 2 2/'(T) + gP-l 

B. Autonomous equations for /(T) dark energy 

To study the dynamics of a general -F(T) model as a dynamical system, 
we introduce the dimensionless variables as follows 

(16) 
(17) 
(18) 

1. (19) 

Using the above relations we can rewrite Eq.0 as the following equation 
Q m = 1 - xi - x 2 - x 3 - x 4 , (20) 



Xl = 


3H 2 ' 


x 2 = 


-2f'(T), 




f(T) 


X3 = 


6H 2 ' 




T 


X4 = 


6H 2 
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with the density parameters 



k 2 pi 

Qi = i = radiation, matter and DE. (21) 

One can rewrite Eqs.(7-ll) as the following equations of motion 

dx\ H N . . 

In = + ( 22 ) 

dxo H 

w = 2mn W< _ < 23) 

W - + < 24 > 

with 

i7 —3x2 — + x i 



# 2 (2m+l)x 2 



(25) 



The autonomous dynamical system Eqs. (22-24) is the general dynamical 
system that describes the cosmological dynamics of f(T) models. 
Where N = ln(-f ) anci 

Tf"(T) 

171 = jrjY) 1 (26) 



From Eq. (|27p one can obtain torsion scalar T as a function of r = then, 
we can say that m is a function of r and m = m(r). 

The effective equation of state Eq. (|15p can be rewritten in terms of x, which 
are defined in Eqs. (16-19) as 

1 -3x 2 - 3x 3 + xi (2m + l)x 2 + 2 , ocA 
" e// = _1 " 3 ( (2m + l)x 2 )( X2 + X3 + I } - (28) 

Eqs. (22-24) show that the results of our analysis depend on X{ Eqs. (16-19) 
and m(r) Eq. (]26p . Therefore, given a form of f(T), we obtain a function 
m{r). For instance, the power-law model f(T) = T + T n corresponds to 
m(r) = n^r- in the (r, m) plane, similar to /(-R) and /(G) [TT]. Also given 



a,i is the initial value of the scale factor. 
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a form of m as a function of r that is cosmologically viable, one may use 
Eqs. ([26|) . ([27|) reconstructing the f(T) viable model. 

On the other hand, from Eqs. (|26p . ()27[) we can derive the following re- 
lation for r: 

d t 

r= r(m + r + 1). (29) 

dN HT K J K ' 

It shows that -^r = if r = and ^y(m + r + 1) does not diverge or 

m + r + 1 = and -§j<r does not diverge. This means that the evolution of 
the system along the curve m(r) stops for any intersection points between 
r = (the m axis) or critical line m = —r — 1 = and the curve m(r). 

3 Phase-space analysis for general f(T) models 

In order to study the dynamical behavior of the system Eqs. (22-24), by 
sitting = we find three continuous lines of critical points with their 
properties in the following relations as: 

Li : (xi = 0, x 2 = -x 3 , x 3 = x 3 ), (30) 
U m = 0, U r = 0, Qde = 1, u e ff = — 1, 



L 2 : (xi =0, x 2 = -2x 3 , £3 = X3, m = 0), (31) 

= X 3 , tt r = 0, Q.DE = 1 — X 3 , UJ e ff = 0, 



L 3 : (xi = x 3 , x 2 = -2x 3 , 2:3 = ^3, = 0), (32) 

VL m = 0, VL r = X 3 , r2z?_E = 1~X 3 , UJ e ff = -. 

It is worth noting that for L±, we have r = —\ and for L 2 and L3 , we 
have r = — 1 corresponding to the definition of Eq. (l27p . Also for L 2 and L 3 
, we have m = so, for these critical lines: 

m (r = -l) = 0. (33) 

Thus for critical lines L 2 and L 3 we have a critical point (r = — 1, m = 0) 
in the (r, m) plane which the system stops in this point in agreement with 
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the critical line m = —r — 1. Note that L2 and L 3 corresponds to scaling 
solutions of the dark energy with matter and radiation respectively. Hence 
in matter or radiation dominated era the m(r) curve should pass through 
the critical point (r = — 1, m = 0) in the (r, m) plane. Also L\ corresponds 
to de-Sitter points (H = 0) which the m(r) curve goes through the critical 
point (r = —5, m = — |). 

The dynamical behavior of f(T) dark energy models have been inves- 
tigated in [22], |23j . Wu et. al, |22j have obtained a critical line for the 
effective dark energy dominated era and two points for matter and radia- 
tion dominated. We will see later that the two points obtained for matter 
and radiation dominated are in agreement with the matter dominated point 
and the radiation dominated point in critical lines L2 and L3, respectively. 
Now we shall consider the properties of any point along the critical lines in 
turn and we study the stability conditions of these points. 

• Pi : de-Sitter point 

(x i = 0, X 2 = -X03, X03 = X03) 

Point Pi is characterized by uj e ff = — 1, ^de = 1 that corresponds to 
de-Sitter solutions (H = 0) for any point along the critical line L\ and 
for any value of £03 • The matrix perturbation eigenvalues for these 
points are 

3 2 lm(m - 1) 8m 3 + 9m 2 + 3m + \ , ^ 

— 4 — ± \ — -A (34) 

2 (2m + l)y (2m + 1) (2m + l) 2 x 03 ' 

with m,} s = m(r = — |) = — |. It is worth noting that the de-Sitter 
point along L\ is stable and an attractor solution when m^ s — > (— ^) _ , 
because the last two eigenvalues are complex and have negative real 
parts. 

Note that, for any value of X03 / Owe have an attractor de sitter point 
with eigenvalues 0, -3, -4 in agreement with [22J. 

• P2 : scaling solutions with matter point 
(xqi = 0, x 02 = -2x 03 , x 03 = x 03 ) 

Point P2 corresponds to scaling solutions with the density parameters 
ratio (jt^: = i_^ 3 03 )• The eigenvalues for these points are given by 

0, -1, 3(l + m, r ), (35) 
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here m jr = The stability of P2 depends on m )7 -. Considering X03 = 
1 this solution leads to P2m = (0, —2, 1) that represents a standard 
matter epoch with Q m = 1, w e // = and the above eigenvalues. 
Hence the necessary condition for P2 to have a standard matter epoch 
is m(r = — 1) = and P2 is a saddle point when 

m r (r = -1) > -1. (36) 

The above relation is the necessary condition to have a standard matter 
epoch which evolve to dark energy dominated era. However, this point 
is stable provided that m r (r = — 1) < —1 . This means that the 
evolution of the system from this point stops and can not move away 
from it and the system does not evolve to dark energy dominated 
phase. 

• P3 : scaling solutions with radiation point 

(xoi = X 3, X 2 = -2x 3, XQ3 = X03) 

At this point, dark energy mimics the evolution of radiation which 
gives the constant ratio = i^f 03 ) an d the point P^ r = (1, —2, 1) 

represents the dominated radiation era with U r = 1 and oj e f f = \- 

0, 1, 4(l+m, r ). (37) 

Note that to have an unstable radiation point m jT .(r = —1) > —1. 

It is worth noting that, P% r and P2m points are in agrement with 
radiation dominated point A and matter dominated point B in [22], 
respectively. They have obtained only two points as radiation and 
matter dominated but, we have obtained two critical lines as scaling 
solutions with radiation point and scaling solutions with matter point. 

We know that a cosmologically viable trajectory would start from the 
radiation era, then pass from a matter dominated epoch followed by 
an accelerated expansion which in the phase space would start from 
P3 then pass the matter dominated point P2 and at last land at the 
stable point P\. In the (r, m) plane, in order to have P3 and P2 the 
curve m(r) should pass the point (r = — 1, m = 0) and the criterion 
for the existence of de-Sitter point Pi is given by passing the curve 
m(r) from point (r = —5, m = —\)- 
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Figure 1: The m(r) curve (red line) for power-law model which satisfies 
the conditions of cosmological viable trajectory, the curve goes through to 
the radiation and matter point (r = — l,m = 0) and the attractor de sitter 
point (r = — i,m = — ■*) and m >r (r = —1) = — ^ > —1. 

4 Specific cases and numerical results 

In this section we study the cosmological viable trajectory for several f(T) 
models by writing m as a function of r and we obtain most of the relevant 
properties of these models. 

I. /(T)=T- F f F 

This model was proposed in |21j-|24j to explain the late-time accelerated 
expansion without including the dark energy component. Also Bengochea 
et. al, [J2] have investigated the observational information for this model by 
using the most recent SN Ia+BAO+CMB data they found that the values 
lie in the ranges n G [0.23,0.03] and £l m 6 [0.25,0.29]. The model presents 
radiation era, matter era and late acceleration phases as the last three phases 
of cosmological evolution in standard model. 

Using Eqs. (j26p . (j27|) one can obtain m(r) = — n(l + -) which is independent 
of a. Note that the necessary condition for the existence of the radiation 
and matter points m(r = —1) = satisfy by this model. m(r = — ^) = — \ 
is the necessary condition to have a stable de sitter point then n = — ^ 
gives a stable de sitter point corresponding to f(T) = T — a^J—T which was 
proposed in |251 [26] explaining the accelerated expansion at late-time. The 
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f(T) 
models 


m(r) 




m(r) = 
r = -1 


m t r > — 1 
r = —1 


m(r) = 
r = -i 

' 2 


T a e bT 

b 

T a er 
r a (ln( 7 T)) fe 


a—r^ 


a+r 2 


a = 1 

a = 1 

11±6V3 
a — 13 


a > -2 
a > 
-1 


a = i 
o= | 

ll±6\/3 
" — 13 


r 

r 2 +2r+a 
r 

1-b 2 , a 2 
b ' "•" fer 

, 2a- 6 
b 


-1 + ^ 



Table 1: Cosmological viability conditions on a number of /(T) models. 
Columns are as follows: 1. f(T) models; 2., 3. Curve m(r) and its first 
derivative m >r corresponds to these models; 4. The necessary condition to 
have radiation and standard matter era; 5. The necessary condition to have 
a standard matter epoch which is followed by the dark energy dominated 
phase, 6. The necessary condition to have a stable de sitter point. 

derivative m with respect to the r is given by 

1 . . 

m >r = n—. (38) 

Hence for this model m jr (r = —1) = — ^ > —1 satisfying the condition 
Eq. ([36|) to have the matter epoch which is followed by the de sitter dark 
energy dominated and is shown in FigfTJ 

II. Exponential model and logarithmic model 

In this section, we consider three cases of the exponential model and 
the logarithmic model which are similar with those of f(R) models and we 
can compare them with the models proposed in f(R) theories [11]. These 
models have been presented in [23]- [26] and we examine the conditions for 
the cosmological viability of them. The numerical results are summarized 
in Table [TJ 

We see that, two exponential models have a matter epoch which is fol- 
lowed by the de sitter point similar to f(R) = Rexp(j^) in f(R) theories [11] 
(see FigJ2]).The phenomenological model with a logarithmic form (m r = — 1) 
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ffiT(Lt}8T) r + TLn(7") 




-10 12 -2-10 



(a) (b) 

ff Tex P m ffT(exp^7) 




-1 (1 1 2 -2-10 




-1 1 2 -2-1 1 2 



(e) (f) 

Figure 2: This figure shows several possible m(r) curves (red lines) for a 
number of f(T) models. Only the case b: f(T) = T + TlnT the m(r) curve 
does not stop the critical line m(r) = —r — 1 at r = — 1 and these models do 
not have the standard matter and the radiation era. In the other cases, the 
m(r) curves intersect the critical line at the matter point (r = — l,m = 0) 
but in all cases the condition Eq. ([36|) does not satisfy to have the standard 
matter era which is followed by the de sitter dark energy dominated and 
system stay in the stable matter point where the values of m, r at r = — 1 
are (—1, —2 and — |) for the cases of a, c and d, respectively which are smaller 
than —1. For the cases e and / the values of m r at r = — 1 are larger than 
— 1 but do not pass through the de sitter point. 
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can not satisfy the condition Eq, (|36p . In this case, unlike the standard 
model, system stays in standard matter era and does not land to the stable 
de sitter point. These models are ruled out by the history of our universe(see 
Fig©. 

5 Conclusion 

The f(T) theory, obtained from generalizing teleparallel gravity, is an al- 
ternative modified gravity to explain the present cosmic acceleration with- 
out including the dark energy. In the present work we investigated the 
dynamical behavior of the recently proposed scenario of f(T) dark en- 
ergy. Performing a detailed phase-space analysis of f(T) dark energy mod- 
els, we obtained three critical lines corresponding to the de sitter point 
(Q m = 0, ri r = 0, VLde = l,^ e ff = —1), scaling solution with matter 
(Q m = l,Q r = 0,Q,de = 0,co e ff = 0) and scaling solution with radi- 
ation (Q m = 0,S7 r = 1,Qde = 0, 0J e ff = 3)- We found conditions for 
the cosmological viability of these models as geometrical constraints on the 
derivatives of these models (r = — T jjp?' , m = ^ JP ). Having de sitter 

point led to passing the curve m(r) into point (r = — \,m = — 5) and to 
have radiation and matter point the curve m(r) should pass through the 
point (r = — 1, m = 0). The stabilities for these points were studied and we 
obtained two kinds of cosmologically viable trajectory which (i) the universe 
would start from an unstable radiation point, then pass a saddle standard 
matter point which would be followed by accelerated expansion de sitter 
point, (ii) The universe started from a saddle radiation epoch, then fell 
onto a stable matter era and the system could not evolve to the dark energy 
dominated epoch. Finally, we investigated the dynamical behaviors for a 
number toy models of f(T) dark energy which were proposed in the more 
literature, recently. These models can exhibit the radiation era, the matter 
era and the late acceleration also, we can compare them with those of f(R) 
theories. 
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